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ABSTRACT 

The recently discovered gamma-ray flares from the Crab nebula are generally at- 
tributed to the magnetic energy release in a highly magnetized region within the 
nebula. I argue that such a region naturally arises in the polar region of the inner neb- 
ula. In pulsar winds, efficient dissipation of the Poynting flux into the plasma energy 
occur only in the equatorial belt where the energy is predominantly transferred by 
alternating fields. At high latitudes, the pulsar wind remains highly magnetized there- 
fore the termination shock in the polar region is weak and the postshock flow remains 
relativistic. I study the structure of this flow and show that the flow at first expands 
and decelerates and then it converges and accelerates. In the converging part of the 
flow, the kink instability triggers the magnetic dissipation. The energy release zone 
occurs at the base of the observed jet. A specific turbulence of relativistically shrinking 
magnetic loops efficiently accelerates particles so that the synchrotron emission in the 
hundreds MeV band, both persistent and flaring, comes from this site. 

Key words: acceleration of particles - pulsars: general - supernova remnants - 
ISM:individual:the Crab Nebula - (magnetohydrodynamics) MHD 



1 INTRODUCTION 

The unexpected discovery of strong, short gamma-ray 
flares from the Crab nebula (|Tavani et al.ll201ll ; lAbdo et al.l 
l201ll ) called attention to the processes that occur in pul- 
sar wind nebulae once again. During these events, the 
flux in a few hundred MeV band grows significantly 
and varies at a time-scale of a half of a day. The ob- 
served properties of the flares place severe limits on pos- 
sible models. The short time-scale implies a very com- 
pact emission region; the high energy of emitted pho- 
tons is an evidence for an extremely efficient accelera- 
tion process. The underlying physical m e chanism is gen- 
erally attributed feednarek fc Idea 2011 ; Uzdenskv et al.l 
201ll;ICerutti et alj|2012t IClausen-Brown fc Lvutikovll20ia 



Sturrock fc Aschwandenl 120121 ') to a rapid magnetic energy 



netic field is kink unstable and would not necessarily pinch 
the flow as much as would otherwise be supp osed; this con- 
jectu re is supported by numerical simulations (|Mizuno et al.l 
l201ll ). In this case the wind magnetization might not be so 
unreasonably small as spherically and axi- symmetric mod- 
els suggest. 

According to the pulsar wind theory, the Poynting 
flux could be converted to the plasma energy on ly via 
dissip ation of variable fields (see, e.g., reviews by lArons! 
(|2007t) and iKirk et al.1 (|2009h >. In the equatorial belt, the 
pulsar magnetic field changes polari ty every half of period 
so th at a striped wind is formed l|Michell Il97ll ; ICoronitil 
1990). The alternating field decays eit her already in the 
wind o r at the terminat io n sho ck jLvubarskv fc Kirkl 



release via, e.g., reconnection, which assumes a magnetically 
dominated region within the nebula. 

Even though the pulsar winds are highly magnetized, 
it is widely believed that at the termination shock, the 
magnetization is already extremely small. This conclu- 
sion is base d on the results of s pherically or axisymmet- 
ric models dRees fc Gunnl 11974 [ Kennel fc Coronitil " 



200 ll; [Kirk fc Skiaeraasenl 120031; Ipyubarskvl l2003b[ 



Petri fc Lvubarskvl 120071 ; IZenitani fc Hoshind 120071 : 



Sironi fc Spitkovskvl l201ll ) therefore in the equatorial 



iKomissarov fc Lvubarskv |2004| ; Del Zanna et al 



il 1 1984 
l2Q0j), 



which demonstrate that the magnetic hoop stress could 
easily distort th e nebula beyond the observational limits. 
iBegelmanl (Il998l ) argued that beyond the shock, the mag- 



belt, a weakly magnetized plasma is injected into the 
nebula; it is this plasma that forms a bright X-ray torus. At 
high latitudes, the magnetic field does not change sign. The 
obliquely rotating magnetosphere excites fast magn etosonic 
waves in this part of the wind (|Bogovalovll200ll ). which 
decay v ia non-linear steep ening and formation of multiple 
shocks (iLvubarskvl l2003af ) . but the fraction of the energy 
transferred by the waves is not large therefore even after 
the waves decay, the flow remains highly magnetized. In 
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this paper, I study the fate of the highly magnetized plasma 
injected into the nebula at high latitudes. 

The paper is organized as follows. In the next section, 
I analyze the lateral distribution of the pulsar wind param- 
eters. In sect. 3, the equations governing the high latitude 
flow in the nebula are derived. In sect. 4, the structure of the 
flow is found. In sect. 5, I shortly discuss stability of the flow 
and mechanisms of the magnetic energy release and particle 
acceleration. Conclusions are presented in sect. 6. 



2 PRELIMINARY CONSIDERATIONS 

The pulsar wind is highly anisotropic; one can conveniently 
adopt the later al distribution of the Poyntin g flux in the 
monopole wind (|Michellll973l ; lBogovalovlll999l l 



F — Fo sin 



(1) 



where F is the energy injected per unit time and unit solid 
angle, 6 the polar angle. Most of the energy is transferred 
in the equatorial belt where the flow has a structure of the 
striped wind. If the pulsar inclination angle is a, the striped 
wind zone is formed at the polar angles exceeding 



(2) 



Since the alternating field decays, a weakly magnetized 
plasma is injected into the nebula in the equatorial belt. 

At high latitudes, 8 < 8o, the magnetic field does 
not change polarity; the obliquely rotating magnetosphere 
just e xcites fast mag netosonic waves propagating out- 



wards (|Bogovalovll200ll '). These waves eventually decay via 
non-linear steepen ing and formation of multiple shocks 
jLvubarskvl l2003aT ) therefore within the nebula, only the 
mean field survives, the energy of variable fields being con- 
verted into the plasma energy. One can calculate the final 
(after the waves decay) wind magnetization assuming that 
the wave has a sine-like form such that the magnetic struc- 
ture is locally presented as 



B = B [l + £sin(fir 



(3) 



where f2 is the pulsar angular velocity, £ < 1 the relative 
amplitude of the wave, Bo the slowly varying (oc 1/r) mean 
field. Throughout the paper, the speed of light is taken to 
be unity. In the relativistic wind, the total Poynting flux is 



F 



(B 2 ) B% 



4tt 



47T 



1 + 



After the waves decay, the Poynting flux decreases to 



F 



(Bf = Bl 

47T 47T ' 



(4) 



(5) 



the residual being converted into the plasma energy. The 
plasma magnetization, defined as the ratio of the Poynting 
to the plasma energy fluxes, after the wave decay is found 
as (one can safely neglect the initial plasma energy) 



F 



F-F 



2_ 



(6) 



One sees that after the waves decay, the wind magnetization 
at high latitudes, 6 < 6o, remains large, a > 2. 



The lateral dependance of the wind magnetization is 
very important beca use it leads to the observed disc-jet di- 
chotomy in PWNe l|Lvubarskvl 120021 ): the disc is formed 
by the weakly magnetized plasma from the striped part of 
the wind whereas the high latitude flow is collimated by 
the magnetic hoop stress i nto a jet-like feature. In simula- 
tions of the PWN structure dKomissarov fc Lvubarsky|2004l ; 
iDel Zanna et al.ll2004 l 20061 ). the magnetization was chosen 
to be zero at the equator and growing with latitude, as it 
should be, but even at high latitudes, a was taken to be well 
below unity. The reason was that high and even moderate 
<t flows are subject to numerical instabilities and moreover 
in axisymmetric simulations, one has to reduce magnetiza- 
tion in order to suppress inappropriately strong elongation 
of the nebula. The study of the PWN structure with re- 
alistic magnetization require s three-dimensional simulations 
so th at the kink instability l|Begelmanlll998l ; iMizuno et aJj 
l201lf) could be taken into account. Recently iKomissarov 
(|2012f ) also pointed out that the residual pulsar wind mag- 
netization is larger than the available simulations assume, 
especially at high latitudes, so that magnetic dissipation via 
three-dimensional effects is crucially important. 

As a preliminary step, I consider a highly magnetized 
conical flow injected into a medium with the constant pres- 
sure po- The pressure within the nebula is determined by 
the energy injected by the equatorial striped wind, which 
transfers most of the energy of the flow. Since the resid- 
ual magnetization of the striped wind zone is small, it is 
terminated at a strong shock. The equatorial radius of the 
termination shock, a, may be estimated from the pressure 
balance condition 



to = p a 



(7) 



Inasmuch as the energy flux in the wind decreases with lat- 
itude, the shock is elongated. In all previous works, the ter- 
mination shock was assumed to be strong everywhere; then 
the distance from the pulsa r to the shock at high latitudes is 
estimated as z = (l/2)9 2 a (jLvubar skv 2002). Now I assume 
that at 6 < Go, the flow is highly magnetized and therefore it 
could be terminated only at a weak shock, which means that 
the postshock flow remains relativistic and radial. Then the 
shock should arise even closer to the pulsar in order to pro- 
vide enough space for the flow to be adjusted to the external 
pressure. 

The postshock flow could be matched to the external 
medium if it is causally connected. In fact the role of the 
termination shock is just to make the flow causally con- 
nected. It is quite possible that a single shock is unable to 
do the job; then a few shocks arise. For example, in weakly 
magnetized models, the high latitude flow passes two shocks: 
a (highly oblique) termination shock and then a rim sh ock 
l|Komissarov fc Lvubarskvll2004 IDel Zanna et al.ll2004T >. In 
any case, the flow eventually becomes causally connected. 
An important point is that for highly magnetized flow, it is 
not necessary to find the shape and number of the shocks 
because the postshock flow anyway remains radial so that 
one can look for the structure of the postshock flow just as- 
suming that it is causally connected and at the inlet of the 
flow, the plasma moves radially. 

Within the nebula, one can neglect the poloidal mag- 
netic field as well as the azimuthal plasma velocity so the 
flow could be considered as composed from magnetic loops. 
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Let the flow with the Lorentz factor 7 subtend the angle 
3> I/7. Then one can consider the flow as being composed 
from magnetic loops moving along the axis and expanding. 
In the frame moving together with the loop upwards with 
the velocity ccos#, the loop expands relativistically with 
the Lorentz factor 7 cxp = #7. The flow remains causally 
connected provided the expansion velocity exceeds the fast 
magnetosonic velocity. In high-<r flows, the last is close to 
the speed of light, the corresponding Lorentz factor being 
7fms = yfa. Then the causality condition, 7 cxp < 7f ms , is 
written as 



6I7 < y/a. 



(8) 



Note that the flow could remain super-magnetospnic and 
even a highly super-magnetosonic if the flow opening angle 
is small. However, if the condition @ is satisfied, the flow 
could be adjuste d to the outer boundary conditions (e.g, 
iBogovalovl (|l997h V Below I assume that the injected radial 
flow already satisfies this condition. 



3 EQUATIONS DESCRIBING THE 
HIGH-LATITUDE FLOW 

Let us for a while consider an axially symmetric flow. One 
can find the flow equations f rom general asymptoti c equa- 
tions for axisymmetric flows l|Lvubarsky||2009l l201ll ). How- 
ever, it is not difficult to derive such equations straightfor- 
wardly taking into account from the beginning that the flow 
within the nebula is purely poloidal whereas the magnetic 
field is purely toroidal so that one can write 



vl, 



B = Be« 



(9) 



Here v is the flow velocity in units of c, 1 the unit vector 
along the flow line, e v is the azimuthal unit vector. Now the 
condition of the flux freezing takes the form 



E = -v x B = -vBn, 



(10) 



where n = 1 x e v is the unit vector in the transverse direc- 
tion. 

In the steady state, one can introduce the electric po- 
tential 



E 



|V*|n. 



(11) 



The plasma flows along the equipotential surfaces therefore 
one can consider f as a specially normalized stream func- 
tion; then the continuity equation is written as 



2nrpyv — 7j(*)|V*|, 



(12) 



where the function r?(^) describes the distribution of the 
mass flux at the inlet of the flow. 

The equation of motion of the hot postshock plasma is 
written in the steady state as 

p 7 (v ■ X7)hjv = -Vp + — [(V ■ E)E + (V x B) x Bl ; (13) 

4-7T 

where p is the plasma proper density, h specific enthalpy. 
One can conveniently make the projections of this equation 
on the direction of the flow, 1, and on the normal to the flow 
lines, n. 

The transfield force-balance equation may be obtained 
by taking the dot product of Equation (|13|) with n. After 



some algebra (see, e.g. jLvubarskv fc Eichlerl |200 ll )). we ob- 
tain in cylindrical coordinates 



1 2 2 E' 

hfnv + — 



where 1Z is the local curvature radius of the flow line (defined 
such that it is positive when the flow line is concave so that 
the flow is collimated towards the axis), 



-n(l- V)l = n-[lx(V x 1)] 



(V x 1). 



(15) 



The flow velocity should be determined from the 
Bernoulli equation representing the energy conservation 
along the flow line. This equation may be obtained by tak- 
ing the dot product of the momentum equation (|13p with 
the longitudinal vector 1. However, one can instead simply 
write down the conserved energy flux per unit mass flux 

hpy 2 v + EB/4 



/*(*); 



which yields, with account of Eqs. (fT0|) . (fTTj) and (JT2J , 
HVtf I 



2r)v 



= Ai(*)- 



(16) 



(17) 



Note that the magnetization parameter a, defined as the 
ratio of the Poynting to the matter energy flux, is presented 
as 



(18) 



Now making use of Eqs. (|10l) . (|lll) and (|17|) one reduces 
the transfield equation (|14|) to the form 



T]pv\V^\ 

2-ktTL 



8irr 2 



7 



1 



+ n • Vp; 



(19) 



In the relativistically hot medium, p oc p 4 ^ s , h = 4p/p. With 
account of the continuity and Bernoulli equations, one can 
write 

4/3 

(20) 



P 

pin 



rfnlh 



H-hrf 



h_ 

hin 



r 2 -y p — fen7ln 
rf n 7i n p - h-y 



r 2 7 p - 



1/3 



(21) 



where the index "in" is referred to the plasma parameters 
at the inlet of the flow. 

Equations JT7]), HU, (f2"U)) and (f2"Tj) form a complete set 
of equations. They should be complemented by appropriate 
boundary conditions. If the flow is confined by the pressure 
of the external medium, po, the pressure balance condition 
should be satisfied at the boundary: 

1 



(B — E ) + p = po. 



(22) 



8tt 

Making use of Eqs. (|10|l and (|11[) . one writes the boundary 
condition in the form 



\ Mi 2 - 1) 



+ p 



PO, 



(23) 



where ^0 is the potential of the last flow line. The boundary 
condition at the axis of the flow is just ty(r = 0) = 0. 

The integrals of motion p and 77 are determined by the 
structure of the pulsar wind. The flow in the wind is radial 
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so that 9, r\ and \i depend only on the polar angle 9. Since 
the wind is highly relativistic, one can present the Poynting 
flux assuming B — E = |V^|, which yields 



F 



4tt 4tt V 99 ) 



(24) 



Here D is the radial distance from the pulsar. Comparing 
this expression with the angular distribution of the Poynting 
flux in the pulsar wind (Q]), one finds the angular distribution 
of the electric potential in the wind: 



* = V47r_Fb(l - cos( 



(25) 



Within the wind, one can write Equation (|17p for the energy 
integral as (note that the wind is cold) 

sin 9 <9* 



fi = 70 + 



2r/ 89 ' 



(26) 



where 70 is the Lorentz factor of the pulsar wind near the 
axis. Eliminating 9 with the aid of Equation (|25|) . one can 
write p(^) close to the axis in the form 



(27) 



which represents in fact the firs t two terms in th e universal 
expansion of in small * (|Lvubarsky||2009l ). The first 

term represents the plasma energy flux whereas the second 
the Poynting flux. Since the pulsar wind is Poynting domi- 
nated, the plasma term could play the role only extremely 
close to the axis. We are interested in the Poynting domi- 
nated domain of the wind, $ S> 7770; therefore we will ne- 
glect the first term Equation (|27|) . 

In the Poynting dominated domain, E 2 /4iyj 2 ^> hp — 
4p, therefore one can also neglect the pressure term in the 
transfield equation (|19|) and in the boundary condition (1231) . 
Since the flow is relativistic, one can also neglect unity as 
compared with j 2 in the right-hand side and take v = 1 in 
the left-hand side of the transfield equation. Now the trans- 
field equation takes the form 



*|V*| 

n 



4r 



-n V 



|V*| 



In the same approximation, Equations ([20 
reduced to 



(28) 



and (1211) are 



h = h m \ — 

,Pin 



1/4 



r?n7in 



1/3 



(29) 



4 THE HIGH-LATITUDE FLOW; 
AXISYMMETRIC CASE 

In principle, the set of equations (|17[) . (|28p and ()29p could be 
reduced to a single equation for 9 just by expressing 7 via 
9 from Equations (|17[) and (|29|) and substituting the result 
into Equation (|28|l . The problem is that in Poynting dom- 
inated flows, 7 is presented in the Bernoulli equation ()17[) 
as a small difference of two large terms, which makes the 
obtained equation for <I> inappropriate for approximate so - 
lution. One can circumvent this difficulty (|Lvubarsky||2009h 
noticing that 7 could be easily found from the transfield 
equation (|28p provided the shape of the magnetic surfaces, 
^(r, 2), is known. An important point is that, in this case, 
an extra accuracy is generally not necessary because in the 



transfield equation, 7 is not presented as a difference of large 
terms. A special care should be taken only if the flow be- 
comes nearly radial because the curvature of the flux sur- 
faces (the left-hand side of Equation (|28p ) is determined 
in this case by small deviations of the flow lines from the 
straight lines. 

Let us for a while neglect corrections of the order hy/fj, 
(i.e. of the order of 1/cr) to the shape of the flux line; validity 
of this approximation will be checked a posteriori. Then the 
Bernoulli Equation (|17p is reduced, with account of Equa- 
tion ((27)) , to 



r|V#| = 24- 



(30) 



which could be considered as an equation for 9. 

One can find an analytical solution to the set of equa- 
tions (|28p and (|30p assuming that the flow subtends a small 
polar angle so that z 3> r. In this case one can take 
n • V = 8/ dr. When looking for the shape of the magnetic 
surfaces, one can conveniently use the unknown function 
r(ty,z) instead of 9(r, z). Then, e.g., 



E= |V*| 



9£ 

8r 



/ 8r \ - 1 



(31) 



In the same approximation, the curvature radius may be pre- 
sented as (note that 1Z is defined to be positive for concave 
surfaces) 



1_ 



8 2 r 
8z 2 ' 



(32) 



Now the transfield equation (|28|) and the Bernoulli equation 
(1301) are reduced to 



8\_ 

8z 2 



8_ 

9* 



9 s 



2* 



8r 



(33) 



(34) 



correspondingly. With account of Equation (|30p the bound- 
ary condition ((23} takes the form 



2np , 



(35) 



The general solution to Equation ((34)) may be presented 



(36) 



where R(z) is an arbitrary function, which is in fact the 
cylindrical radius of the flow. In order to find the function 
R(z), let us substitute Eq. ((36p into the left-hand side of Eq. 
and integrate the obtained equation from to 1 3E r c£El : 



1 Formally speaking, one could not integrate from since the 
solution (I36D was obtained only in the Poynting dominated do- 
main, where the first term in the full Bernoulli equation l|17l) is 
neglected. Therefore this solution becomes invalid close enough 
to the axis where the flow ceases to be Poynting dominated 
(the Poynting flux goes to zero at the axis, see Equation ( 1271 ). 
However, the integrands grow with ^ so that the lower bound- 
ary condition cou ld be continued to zero. For more details, see 
llLvubarskvll200gF) where a solution for the \1/ <g ^0 part of the 
flow, including the matter dominated core of the jet, has been 
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3 dz 2 



(7 2 )*=>i< (J 



(37) 



Making use of the boundary condition (135[l one gets a simple 
equation for R: 



where 

67TJJ0 



P = 



(38) 



(39) 



Equation (J38J is a partial case of the governing equation 
that determines a ll the properties of highly magnetized flows 
|Lvubarskvll2009h . 

Solution to Equation (|38|l describes the postshock flow 
at high latitudes where the wind remains highly magnetized. 
Let the high magnetization flow subtend a small polar angle 
80 ■ Then one can write, with account of Equations Q and 



P = 



6 



(40) 



Since the shock is weak, the flow just beyond the shock 
remains nearly radial so that if the flow line = \]/o enters 
the shock at the distance zo from the pulsar, the solution to 
Equation (|38p should satisfy the conditions: R(zo) = zo8q; 
R'(zq) = 8q. Such a solution is simply 



R = — -= sin y pz = — = 



9la . y/Ez 



V6 



(41) 



Note that the solution is independent of the position of the 
shock, zo. This is because the postshock flow remains any- 
way radial. At the end of this section, I present an estimate 
for z Q . 

The Lorentz factor of the flow is obtained by substi- 
tuting the solution (|36f) and (|41[) into Equation (|33[) and 
integrating from zero to ^; this yields 



7 = 



3^o 



la sm 



V6z 
9 2 a 



(42) 



One sees that initially the flow expands and decelerates and 
then converges and accelerates. Substituting the obtained 
solutions into Equation (|29p . one finds 



h oc sin 



-1/3 



(43) 



which means that the flow is cooled when it expands and 
heated again when it converges. 

The flow is maximally expanded at the distance 



zi = fy a. 



2^6 



(44) 



At this distance, the radius and the Lorentz factor of the 
flow are 



Ri = —0 zi 

7T 



V6 



7i 



v/3 



(45) 



obtained; this solution is matched smoothly with the solution in 
the Poynting dominated domain. 



Note that if the flow were weakly magnetized so that the 
termination shock were stro ng, the shock wo uld arise at the 
distance z = (l/2)#oa ~ Zi 

l|Lvubarskvll2002l L roughly equal 
to z\ . In a highly magnetized flow, the shock is much closer 
to the pulsar. 

The above solution was obtained by neglecting the hj 
term in the Bernoulli Equation (|17[) . This means that the ex- 
pression ()36[) for the shape of the flow lines, where R(z) sat- 
isfies the equation (|38|) . is valid to within /ry/^ = 1/a <C 1. 
This approximation fails if the flow lines are nearly straight 
because in this case, the curvature of the flow lines (the left- 
hand side of the transfield Equation (|33p ) is determined by 
small deviations from the straight line so that even small 
corrections to the flow line shape could not be neglected. In 
our solution, this happens near the origin of the flow, z ~ 0, 
and near the converging point, z « 2z\. 

In order to find limits of applicability of the solution, 
let us present the shape of the flow lines as (cf. Equation 



R(r) 



(1 + 5), 



(46) 



where <5(>t,.z) describes a small correction to the shape of 
the flux surfaces due to a nonzero fvy/fj,. Substituting this 
expression into the full Bernoulli Equation (|17[) and expand- 
ing, with account of Equation (|27[) and (|31[) . in small S, one 
gets 



(47) 



2<S> 2 05 

which implies S ~ h'y/fi = 1/a. One can neglect this correc- 
tion provided the contribution to the curvature of the flow 
line due to this correction, 



d 2 RS 



dz 



RS 
75" 



R 



(48) 



remains less than d 2 R/dz 2 . According to equation (I38|l . the 
last is just equal to —/3R therefore the condition that one 
can neglect the corrections of the order of 1/a to the shape 
of the flow line, which is in fact the validity condition for the 
solution described by Equations (|36|l and (|41|l . is written as 



z^>—= = —^a= — — Zl . (49) 



Taking into account that according to Equation (1421) . the 
Lorentz factor of the flow at small z is presented as 7 = 
9oa/(i/2 z), one sees that this condition may be written as 
0o7 "C y/a, which is in fact the condition (j8} of the causal 
connection of the flow. Therefore our solution is valid only 
in the domain of the causal connection of the flow. 

The full structure of the flow could be formally ob- 
tained by matching the above solution with the free pul- 
sar wind by inserting shocks and making use of the shock 
jump conditions. I do not address here the entire problem 
but just assume that shock(s) make the flow marginally 
causally connected, after which my solution is valid. This 
conjecture could be justified by contradiction: if the flow re- 
mains causally disconnected (even having passed a shock), it 
could not be adjusted to the outer conditions and therefore 
more shocks must arise; if the flow is causally connected, it 
is smoothly adjusted to the outer boundary conditions ac- 
cording to the above solution. Therefore one could expect 
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that the shock(s) just make the flow marginally connected, 
after which the above solution works. This means, according 
to the estimate (|49jl . that the shock(s) arise approximately 
at the distance 

zo = —;=a « -=. (50) 

According to this solution, the flow is focused to a point 
2 = 2zi at the axis. Close enough to this point, just like 
close to the origin, the solution becomes invalid because the 
flow lines become nearly straight. But what seems to be 
more important is that such a converging flow ceases to be 
axisymmetric because of development of instabilities. This 
issue is discussed in the next section. 



5 INSTABILITIES, TURBULENCE AND 
PARTICLE ACCELERATION 

The axisymmetric flow considered in the previous section 
could not remain axisymmetric. The flow is composed from 
magnetic loops disconnected one from another therefore 
they could easily come apart. In the volume of the flow, 
the axisymmetry could be destroyed by the kink instabil- 
ity |Begelmanlll998l ; iMizuno et ail 1201 ll ~). Near t he bound- 
ary o f the flow, the Kelvin-Helmholtz instability (Begclman 
I1999T ) develops. In relativistic flows, these instabilities de- 
velop slowly because of relativistic time delay. However, if 
the flow converges, even small perturbations eventually de- 
stroy the regular structure. If two converging loops initially 
shifted one with respect to another by a displacement much 
less than their radius, the distortion becomes strong when 
the radius approaches the initial displacement. Therefore 
when the axixymmetric flow is focused into a point at the 
axis, the magnetic loops come apart close enough to the 
converging point so that a specific turbulence of shrinking 
magnetic loops emerges. 

Taking into account that the hoop stress within the loop 
is not counterbalanced by either the poloidal magnetic field 
or the plasma pressure, the loops shrink with relativistic 
velocity, 7 ~ y/a, until the plasma energy reaches the mag- 
netic energy. Therefore one concludes that the energy of the 
high latitude, highly magnetized flow is released close to the 

converging point , z 2zi. In this picture, the observed jet 

I Weisskopf et aUl200Ch begins in the vicinity of this point. 
One can expect also an efficient particle acceleration be- 
cause in such a turbulence, the electric and magnetic fields 
are nearly equal. 

The observed synchrotron spectrum from the Crab is 
extended up to ~ 100 MeV, which implies particle accel- 
eration on a time scale of the Larmor gyration period so 
that the acceler ating electric field should be as strong as th e 
magnetic field (|Guilbert et all 1 19831 ; Ide Jager et aiTll99fJ) . 



The recent discovery of da ylong gamma ray flares in the 
band of a few hundred MeV jT avani et al.ll201ll;lAbdo et alJ 



2011 
2012 


;IStriani et al]|201ll;IVittorini et al.ll201ll;lBuehler et al. 

) poses even a larger challenge to the acceleration mod- 


els ( 


Bednarek & Ided 20 111; Komissarov & Lvutikov 


2011; 


Uzdenskv et al.l2011; Bvkov et al.ll2012l; Cerutti et al.l 


2012; 


Clausen-Brown & Lvutikovl 120121: Isturrock & Aschwanden 


20121). The relativistic turbulence could in principle resolve 



the problem. On the one hand, at relativistic turbulent ve- 



locities or, which is the same, at E « B, even the second 
order Fermi process is efficient enough to permanently accel- 
erate electrons up to the energies sufficient to emit 100 MeV 
synchrotron photons. On the other hand, random relativis- 
tic bulk motions could occasionally produce flares even in 
a harder energy ban d when a local rad iating "knot" moves 
toward the observer l|Yuan et al.ll201lh . 

The efficient particle acceleration could occur not 
only via Fermi mechanism but also via reconnection 
ijUzdenskv et al.ll201ll ; ICerutti et ai1l2012t ). When different 
magnetic loops slide one over another, strong field gradients, 
and therefore strong currents, arise. The magnetic reconnec- 
tion comes into play when the current velocity approaches 
the speed of light. The corresponding scale, S ~ B/(8iren), 
could be expressed via the plasma multiplicity defined as the 
ratio of the plasma density to the Goldreich- Julian density 
in the pulsar magnetosphere, k — encP/B, where P is the 
pulsar period. Taking into account that within the magneto- 
sphere, n oc B, whereas beyond the light cylinder, n oc B/r, 
one finds 5 ~ r/n. This scale may be considered as the dis- 
sipation scale of the turbulence; at this scale the magnetic 
energy is converted into the plasma energy. 

The detailed analysis of the particle acceleration is be- 
yond the scope of this paper; the above consideration just 
shows that a plausible site for the synchrotron emission in 
the hundreds MeV band, and in particular for gamma-ray 
flares, is the region close to the point z = 2zi toward which 
the high latitude flow is focused. According to Equation (|44p . 
this region is at the axis of the system at the distance a 
few times less than the equatorial radius of the termination 
shock. In this region, the energy transferred by the strongly 
magnetized high latitude flow is released. T his region may be 
ident ified with the base of the observed jet (Weisskop f et al.l 
|2000D . 



6 CONCLUSIONS 

In this paper, I considered the fate of the high latitude flow 
in the pulsar wind. This part of t he wi nd remains strongly 
magnetized (see also iKomissaro 3 <|2012l )) therefore it is ter- 
minated at a weak shock (or at a system of weak shocks), be- 
yond which the flow is still radial and relativistic. The shock 
arises very close to the pulsar, much closer than if it were 
strong; the higher the magnetization of the flow, the closer 
the shock to the pulsar. Beyond the shock, the flow still ex- 
pands but decelerates and eventually becomes to converge 
because the magnetic hoop stress is not counterbalance ei- 
ther by the poloidal field or by the plasma pressure. 

In the converging flow, magnetic energy is converted 
into the plasma energy therefore the plasma accelerates and 
heats. Even small perturbations due to, e.g., kink instabil- 
ity eventually destroy the converging flow so that magnetic 
loops come apart and then shrink independently of each 
other producing relativistic turbulence. Hence one can ex- 
pect that the whole energy of the highly magnetized part of 
the pulsar wind is released in a very small region close to the 
converging point, which occurs on the axis of the system at 
the distance from the pulsar ~ O^a, where 8q is the opening 
angle of the highly magnetized part of the wind, a the equa- 
torial radius of the termination shock. I identify this region 
with the base of the observed jet. Relativistic turbulent mo- 
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tions in highly magnetized plasma imply E « B so that in 
the energy release region, particles could be efficiently ac- 
celerated either via the second order Fermi mechanism or 
via the magnetic reconnection. Therefore the synchrotron 
gamma-ray emission in the hundreds MeV band, both per- 
sistent and flaring, could come from a small region at the 
base of the jet. 
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